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We studied the intensity of resonant Raman scattering due to optical phonons in a planar II- 
VI type semiconductor microcavity in the regime of strong coupling between light and matter. 
Two different sets of independent experiments were performed at near outgoing resonance with the 
middle polariton (MP) branch of the cavity. In the first, the Stokes-shifted photons were kept at 
exact resonance with the MP, varying the photonic or excitonic character of the polariton. In the 
second, only the incoming light wavelength was varied, and the resonant profile of the inelastic 
scattered intensity was studied when the system was tuned out of the resonant condition. Taking 
some matrix elements as free parameters, both independent experiments are quantitatively described 
by a model which incorporates lifetime effects in both excitons and photons, and the coupling of 
the cavity photons to the electron-hole continuum. The model is solved using a Green's function 
approach which treats the exciton-photon coupling nonperturbatively. 

PACS numbers: 78.30.Fs, 78.67.-n, 71.36.+C 

Keywords: Raman spectroscopy, II- VI semiconductors, optical microcavities, optical confinement, cavity 
exciton-polaritons 



I. INTRODUCTION AND MOTIVATION 

Semiconductor microcavities (MC's) are micron-scaled 
photonic structures, consisting of two planar distributed 
Bragg reflectors (DBR's) separated by a spacer with a 
thickness of the order of the photon wavelength, forming 
a high quality Fabry-Perot resonator, and leading to a 
localized photonic mode inside the MCii This spacer 
usually embeds quantum wells (QW's) conveniently 
located at the antinodes of the optical-cavities electric 
field. The strong vertical confinement of both the 
electric field and the QW excitons result in a strong and 
by design controllable light-matter interaction. When 
this MC mode is resonant with the excitonic transitions, 
and when the life time broadening of the MC-photons 
(7/) and QW-excitons (/y e ) are smaller than the 
photon-exciton interaction (O), the system is under the 
strong-coupling regime, where coherent Rabi oscillations 
between the bare states occur. This regime is evidenced 
by the characteristic anticrossing behavior exhibited by 
the exciton and the cavity mode when they are brought 
into resonance. As a consequence, the "new" eigenstates 
of the system are strongly mixed exciton-photon states, 
cavity-polaritonsX^^ These quasi-particles have shown 
to exhibit a large number of interesting properties, most 
of them arising from their inherited bosonic character. 
This character together with their very small effective 
mass (typically 10 -4 times smaller than that of the bare 
excitons) enables processes like the condensation into a 
quantum degenerate Bose-Einstein stated in thermal 



equilibrium with the lattice, 7 and the boser effect 
or stimulated emission^^^ among other non-linear 
phenomenal This made the cavity-polariton systems 
not only very interesting from the basic fundamental 
point of view, but also very promising in the field of 
technological application in coherent light emitting 
devices i^ii i 11 : 12 Efforts have been devoted incessantly 
to the study of cavity-polaritons experimentally and 
theoretically during the last fifteen years^^ since the 
first experimental observation* 1 ^ 

The process of inelastic light scattering due to vi- 
brations in semiconductors and in microcavities for the 
case where the energy of the light mode is far from 
the corresponding one of electronic excitations, is usu- 
ally very well described by treating the light-matter cou- 
pling perturbatively 14 i 15 In this case the Raman effect 
is mediated by electronic or excitonic states by a sec- 
ond order process in exciton-photon interaction. In- 
stead, in microcavities when the cavity mode is tuned 
near the QW-exciton energy the perturbative approach 
failsi 16 i 17 ' 18 i 19 ' 20 The exact treatment of the light-matter 
interaction, i.e. the cavity exciton-polariton picture, is 
essential and the inelastic process (now mediated by po- 
laritons) must be described by a first order perturbation 
theory of the polariton-phonon interactioni 16 : 18 ? 21 

Since the Raman process due to optical phonons is a 
coherent process, it is very sensitive to the damping of 
the resonant state j 16 i 17 ' 18 > 19 Resonant Raman scattering 
(RRS) experiments as a function of photon-exciton de- 
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tuning evidence the dependence of the polariton-phonon 
and polariton-external-photon coupling with the mode 
mixing in the process. This is distinctly demonstrated by 
the maximum in the scattering efficiency near the mode 
anticrossing (zero detuning) and the decreasing efficiency 
when the polariton character approaches either the pure 
exciton or pure photon stateJ^- 9 Although existing sim- 
plified theoretical models describe the phenomenological 
and qualitative behavior of the experimental results, the 
relevance of rigorously including polariton dcphasing and 
dumping effect has been strongly emphasizedj 17 i 19 ' 22 

In this paper we present detailed RRS experimental 
results together with a theoretical improvement of the 
existing models that rigorously takes polariton lifetimes 
into account. The experiments are performed on II- VI 
semiconductor based microcavities and consist essen- 
tially of two types: the first showing the behavior of the 
Raman intensity as a function of cavity-photon-exciton 
detuning at exact resonance, and the second showing 
the behavior of the outgoing scattered process for a fixed 
detuning, and performing a scan with the excitation 
energy. Experiments are compared to the developed 
theory that uses a Green's function approach to include 
the damping effects and the electron-hole (free-exciton) 
states in a formal way. Using the matrix elements in 
the theory as free parameters, the theory describes 
accurately the Raman intensity simultaneously for both 
types of experiments. Bare life time broadening are 
obtained from a fitting procedure and the effects of 
the exciton lifetime in relation to the cavity-polariton 
lifetime are discussed. An important improvement to 
the previously presented theoretical approache s 23 ! 24 to 
describe the experimental observations quantitatively is 
obtained in the present work by accounting also for the 
second type of experiment, which provides additional se- 
vere constrains to the unknown parameters of the theory. 

The paper is organized as follows: In section [II] the 
samples are shortly described and the experiments are 
presented. In section IIIII the Experimental results are 
shown. The theoretical model is introduced in section lTVl 
a short description of the existing phenomenological ap- 
proaches is given, and afterwards the model Hamiltonian 
and Green's function approach for the Raman intensity 
are explained. In section [V] the experimental and cal- 
culated results are analyzed and discussed. Finally the 
conclusions of this work are summarized. 



II. EXPERIMENTAL DETAILS 

Experiments are all performed on a II- VI semicon- 
ductor microcavity. CdTe-based cavities are interest- 
ing since excitonic binding energies in CdTe quantum 
wells are much larger (~ 25 meV) as compared to GaAs 
(^10 meV). This is a consequence of a stronger electron- 
hole attraction that makes the excitonic states in this 
compounds much more stable. In addition, the ionic 



character of these compounds leads to a stronger light- 
matter coupling and hence to much larger Rabi gaps and 
stronger polaritonic effects compared to the GaAs-based 
structures. On the other hand Raman efficiency due to 
optical phonons is also larger allowing clearer studies, 
specially when the Stokes-shifted photon is tuned at near 
resonance to a polaritonic state, where luminescence hin- 
ders its observation. 

The studied cavity was grown (lOO)-oriented by molec- 
ular beam epitaxy, and consists of two distributed Bragg 
reflectors (DBRs) enclosing a A/2 Cdo.4Mgo.6Te spacer. 
Three centered CdTe quantum wells of nominal width 
72 A and separated by 69 A from each other, are em- 
bedded in the spacer. The cavity mirrors consisted of 
15.5 periods of Cdo.4Mgo.5Te / Cdo.75Mgo.25Te A/4 lay- 
ers for the top DBR, and 21 pairs for the DBR at the 
bottom, to balance the effect of the air-substrate asym- 
metry. The layers were grown in such a way that their 
thickness decreases in one direction (like in a wcgde) , and 
the cavity mode can thus be tuned by displacing the po- 
sition of the laser spot on the sample. Experiments were 
performed at back-scattering z(x,x')z' geometry, where 
z(z') and x(x') describe the incident (scattered) beam di- 
rection and light polarization, respectively* 2 ^ Here, x and 
x' correspond to the (110) direction, z to the (001) epi- 
taxial growth direction, and since the configuration is 
back-scattering, z' = z. A tunable cw Tksapphire laser 
was used as excitation source at near normal incidence, 
with power < 50 fxW, and focalized to a spot of ~ 50 
fim. The fixed collection cone normal to the sample sur- 
face was ~ 3°, corresponding to an uncertainty in the 
in-plane wave vector of k\\ « 4.2 x 10 3 cm . A liquid- 
He cryostat was used to set the temperature at 4.5 K. At 
this temperature, the sample shows a double anticrossing 
of the cavity mode with both the ElHl(ls) (-Xi s ) and 
the ElHl(2s) (^2s) heavy-hole exciton modes. The dis- 
persion of the middle and lower polariton branches (MP 
and LP, respectively), derived from photoluminescence 
experiments, is displayed in Fig. [TJl as function of the 
spot position x (in mm) on the sample. Unfortunately 
no direct observation of the presence of the upper po- 
lariton (UP) branch was found with luminescence. Using 
a simple three-coupled-mode model, the best fits of the 
experimental MP and LP branches was obtained for a 
Rabi-splitting of Q/_i s ~ 13meV and 0/_2s — 3.5 meV 
between the cavity mode and the Is and 2s excitons, re- 
spectively. The dashed and dash-dotted lines correspond 
to the bare (non- interacting) exciton (Xi s and Al2s)> 
and cavity-photon (/) dispersions used in the model, 
respectively^ The vertical dotted lines indicate the po- 
sitions of the anticrossings, i.e. the minimal distance be- 
tween the MP-LP branches and the UP-MP branches. 
The slight positive slope of the bare exciton dispersions 
is due to the small thickness gradient, that also affects 
the QW's. 

Two different and independent Raman experiments 
were performed observing the inelastic scattered photons 
due to the Raman allowed CdTe QW longitudinal op- 
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FIG. 1: (Color online) (a) Dispersions of the cavity-polariton 
modes as functions of the spot position on the sample, derived 
from photoluminescence. Experimental data: lower (LP, 
empty dots) and middle polariton (MP, solid dots) branches. 
The full lines are a fit to the experiments with a three-coupled- 
mode model, using the plotted bare cavity-photon (/, dashed- 
dotted) and exciton (X\ s and X 2a , dashed) curves, and Rabi- 
splitting of Sl/_i s ~ 13meV and fi/_i s ~ 3.5 meV, as input, 
(b) Photon (|^4/| 2 , full line) and exciton ( | ^4i s | 2 and \A 2a \ 2 , 
dashed and dash-dotted lines) strengths of the MP, derived 



from the fit in panel (a). A ls 



A ls A 2s + A* ls A 2 s (short- 



dashes), is the exciton crossed interference term. 



tical (LO) phonons. These phonons have an energy of 
~ 172 cm -1 (~ 21.3 meV). The first experiment, con- 
sisted in setting the Stokes-shifted photons always at ex- 
act outgoing resonance with the MP, and changing the 
cavity-photon-exciton detuning by displac ing the laser 
spot on the sample, as described in Refs. flllUs ]. The sec- 
ond experiment consisted in keeping the detuning fixed 
(i.e. at a fixed spot position), and varying the laser en- 
ergy at near outgoing resonance, in order to scan the MP 
with the Stokes-shifted photons. In both cases the spec- 
tra were analyzed using a Jobin-Yvon T64000 triple spec- 
trometer in substractive mode, equipped with a liquid- 
s-cooled charged-coupled device (CCD). 



III. EXPERIMENTAL RESULTS 

Typical spectra for the RRS experiment (of the first 
kind described in section [TT|, are shown in Fig. O), for 
different MP mode energies (i.e., laser spot positions). 
The Raman peak is highlighted on top of the MP lumi- 
nescence. For each spectrum the exact resonance was set 
tuning the laser energy in order to maximize the Raman 
peaks intensity. 

In the bottom panel (Fig. the Raman intensity, 
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FIG. 2: (Color online) (a) Experimental and fitted differences 
between the energies of the MP and LP branches (left axis), 
and UP and MP (right axis), as function of the MP energy 
derived form Fig. (b) Raman spectra at exact outgoing 
resonance with the MP for varying detuning. The Raman 
peaks are indicated by the filled areas. The arrows labeled 
as (a) and (b) indicate the equivalent positions of the Raman 
scans of Fig. [3^. and [3)3 respectively, (c) Resonant Raman in- 
tensity (full dots) derived from panel (b). The dashed curve 
corresponds to the calculation using the simple model with- 
out damping effects, and the full curve corresponds to the fit 
using the theory derived in the text. The vertical dotted lines 
indicate the position of both anticrossings. And the horizon- 
tal arrows in panel (c) indicate the shift of the RRS intensity 
maxima. 



derived from the spectra in panel (b), is shown as func- 
tion of the MP energy. As additional information in the 
top panel (Fig. the difference between the measured 
MP and LP energies is also displayed. The full (dashed) 
line corresponds to the calculated differences Emp — Elp 
{Ejjp — Emp), using the fitted curves and parameters 
of Fig. [Iji obtained from the three-coupled-mode model. 
The vertical dotted lines indicate the anticrossing posi- 
tions, i.e. the minima of the respective energy differences. 

In Fig. ^Bp, for low MP energies, luminescence corre- 
sponding to localized excitonic states (labeled as Xi oc ) 
can be distinguished. 18 And for increasing MP energies, 
i.e. for increasing photon- Xi s detuning, the overwhelm- 
ing LP luminescence appears from the left. One thing 
to notice is how the MP luminescence increases when its 
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FIG. 3: (Color online) Outgoing Raman scans for a fixed 
photon-exciton detuning, and varying laser excitation energy. 
Scan (a) corresponds to a spot position of x ~ 2.5 mm and 
(b) to x ~ 3.1mm. The Raman peak moves on top of the 
broad MP luminescence peak and resonates at its position, 
indicated by the vertical dotted line. The resonant situation 
is highlighted with a thicker curve, and corresponds to the 
positions indicated with vertical arrows in Fig. [5] and labeled 
with (a) and (b) respectively. Note that the broad lumines- 
cence peak in (b) is ~ 3 times larger than in (a). 



energy approaches the X2 S exciton at ~ 1.665 eV. Note 
that the luminescence maximum is reached slightly be- 
low this energy. This typical behavior of luminescence 
in the presence of an anticrossing is refered to as "cavity 
pulling", in this case between the MP and the UPi2£ A 
similar behavior is observed for the luminescence coming 
from the LP branch, slightly before the anticrossing of 
the MC photon and the X\ s exciton. 

As for the second type of experiment, typical first order 
LO-phonon Raman scans for varying laser energy close to 
outgoing resonance with the MP at two different spot po- 
sitions (i.e., detunings) are shown in Fig. [31 In Fig. [3^, the 
MP is centered at ~ 1.6548 eV (corresponding to a spot 
position on the sample of about x ~ 2.5 mm), whereas in 
Fig. [3)3 it is centered at ~ 1.6599 eV (for x ~ 3.1mm). 
The spectra corresponding to the outgoing MP resonant 
condition are displayed with thicker lines, and the corre- 
sponding equivalent resonant spectra are indicated with 
arrows in Fig. [2)3. Note that the photoluminescence back- 
ground of both scans in Fig. [3J differ considerably, as is 
clearly observed in the corresponding resonant spectra 
on Fig. [2)3. Figure [4] shows the Raman intensity profiles 
(a) and (b) derived respectively from Fig. [3^, and [3)3. A 
third profile (c) is also displayed, which corresponds to 
an energy of the MP at ~ 1.6533 eV (for a spot position 
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FIG. 4: (Color online) Raman intensity scans for varying laser 
energy close to the MP outgoing resonance for three different 
detuning. The MP energy corresponds respectively to (a) 
E M p ~ 1.6548 eV, (b) ~ 1.6599 eV, and (c) ~ 1.6533 eV. 
And the positions on the sample are (a) x ~ 2.5 mm, (b) 
~ 3.1mm, and (c) ~ 2.2 mm. The Raman intensity resonates 
at the MP energy. The experimental intensity labeled as (a) 
and (b) are derived from the spectra shown in Fig. [3ji and 
[3]b respectively. The full lines are a fit to the data using the 
theory derived in the text (Sec. IIV| l. 



x ~ 2.2mm)p. Note the asymmetry of the intensity 
scans labeled as (a) and (c) in Fig. [4] towards lower en- 
ergies. This might be ascribed to the strong presence of 
the LP. At lower energies the resonant contribution due 
to the LP starts to contribute to the Raman peaks inten- 
sity, thus inducing the asymmetry in the Raman scan. 
Scans (a) and (c) are taken near the anticrossing, where 
both branches (MP and LP) are closer to each other. 
This is not the case for the scan (b) (see Fig. [3Jd and 
[4]). This scan is far from the LP and at this detuning its 
resonant contribution is less strong. Thus this intensity 
scan is more symmetric. 

In the next section the theoretical background will be 
addressed, the Green's function formalism presented, and 
the calculation of the Raman intensity within this model 
will be described. 



IV. THEORETICAL DESCRIPTION 
A. The model 

In a cavity polariton system, only MC confined pho- 
tons and QW confine excitons with the same in-plane 
wave vector (fcii) are coupled^ Since the above described 
experiments are performed in near-backscattering geom- 
etry, i.e. z(x', x')z, fcii is considered to be zero (fen = 0)^ 
In what follows the total wave vector and polarization in- 
dices will be suppressed, in order to simplify the notation. 

The following Hamiltonian, for a given k\\ , is considered 
to describe the cavity exciton-polariton system^ 



H = H f + H e + H c + H a 



(1) 
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The first term corresponds to the photonic part of the 
Hamiltonian: 

H f = E f pf + £ e p r\r v + £(V>t/ + H.c). (2) 

p p 

Here f< creates a MC photon, which is hybridized to 
a continuum of radiative modes described by creation 
operators r' p . The effect of the latter is to provide a finite 
life time (7/) to the cavity mode, even in the absence of 
light- matter interactions. 

Similarly, the excitonic part of H is described as 

H e = E i e Ui + e iQ d l q diq + J2( Vi 1 d i ei + H - C ') ' 

i iq iq 

(3) 

where ej creates exciton states that will couple strongly 
to the cavity photon mode. For example e\ creates a Is 
exciton state of energy E\ s . Analogous to the case of the 
MC photon, we assume that each exciton mode mixes 
with a continuum of bosonic excitations (described by 
the operators d ig ), which broadens its spectral density. 
The effect of this mixing is to introduce a broadening "f e i 
for the excitons even without interaction with the light. 
The detailed description of this broadening is beyond the 
scope of this work. Numerous physical mechanisms are 
known to contribute to the damping of the excitonic wave 
function and hereby to the line broadening. Some of 
them are for example: impurities, defect-scattering, in- 
terface roughness, well- width fluctuations, and alloy dis- 
order, which have been demonstrated to be present in this 
kind of II- VI based microcavities? 2 ^ Whittaker has shown 
that the inhomogeneous broadening introduced by disor- 
der can be described accurately introducing an imaginary 
part in the self energy of the exciton Green function i 29 i 30 
As we show below, H e has a similar effect within our the- 
ory. Therefore the broadenings 7^ [see Eq. ([15])] can be 
regarded as due to lifetime broadening. 

The third term in Eq. (fTJ) describes the electronic ex- 
citations that do not couple strongly to the MC photon 
mode, 

H c = ^ e k ala k ■ (4) 

k 

These states include discrete exciton states as well as 
states of the exciton-continuum (i.e. the continuum of 
electron- hole excitations). The latter correspond to the 
eigenstates of the exciton hydrogen-like Hamiltonian that 
are not bound and have positive internal energy. 31 ' 32 The 
operators, which describe this exciton continuum, can be 
well approximated by a' k — ctufc, where Vk destroys a 
valence electron with 2D wave vector k and polarization 
<7, and c\. creates an electron in the conduction band with 
the same k and a (so that the total wave vector kn = 0). 

Note that the commutator [afe,ol] = vtvk — c \ c k- For 
the temperatures and laser intensities of the experiments, 
all valence (conduction) states are practically occupied 



(unoccupied) and therefore with high accuracy vlvk = 1 
(c\ck = 0). Therefore the operators entering Eq. (|4]) can 
be regarded as ordinary bosons i 33 ' 34 

The separation of the bound excitonic states into those 
included in H e or H c is in principle arbitrary. In particu- 
lar, when the width of the weakly bound excitonic states 
near the bottom of the gap is larger than the separa- 
tion in energy between two next states, these states can 
be well described as a continuum, beginning at some en- 
ergy E' x , slightly below the energy of the first excitonic 
state not included in H e . In our specific treatment for 
the experimentally studied sample, we describe all states 
included in H c as a continuum beginning at E' x , and in- 
clude the two lowest excitonic states (E\ s and i?2s) that 
couple "strongly" to the MC photon in H e . For simplic- 
ity we take E' x = -E^s- 

Finally, the last term of the Hamiltonian Eq. ([I]) de- 
scribes the mixing between light and matter 

Hnxix - + H - c -) + + H.c). (5) 

i k 

According to the quantum theory of radiation, in the 
dipolar approximation Vk is proportional to {k c \p a \k v ), 
where p a is the momentum operator in the direction of 
the electric field, and \k c ) (\k v )) is the wave function 
which corresponds to the operator Ck (vfc)- For simplic- 
ity, we are neglecting terms which do not conserve the 
number of bosons. They can be included without much 
complications, but their effect is small for the energies 
of interest^ 

The Hamiltonian Eq. ([1]) is quadratic and can be for- 
mally diagonalized by a Bogoliubov transformation to 
the form 

H = Y J V vV lv v , (6) 

V 

where the generalized polariton operators p' v , are related 
to the creation operators entering Eq. |T]) by the linear 
combination 

pt=£^/3j. (7) 

3 

(3j stands for any operator entering Eq. (TT]), and A V j are 
the coefficients that give the (detuning dependent) cor- 
responding photonic or excitonic weight of the polariton 
state v. 

Note that if only the first terms in Eqs. ((2|), ([3]) and 
([5]) are kept, the simple standard coupled- mode model 
is obtained. 19 ' 22 For instance, the three-coupled-mode 
model results, when the MC mode couples strongly only 
to the i — Is and the i — 2s exciton states, and the 
values of the interaction are given by 2V\ = f2/_i s and 
2V2 = £l/-2s, as it is the case for the above described 
experiments. 
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B. The Raman intensity 

It has been shown previously ; 18 ! 19 ' 22 ! 23 ' 24 that the out- 
going Raman cross section mediated by cavity polaritons 
can be calculated by first order perturbation theory as 



J(w) oc Ti W l ^ s T s 



(8) 



where Tj describes the incoming channel, i.e. the con- 
version of the incident external photon to the initial po- 
lariton state inside the sample. For the case of outgoing 
resonance and the scattering geometry used in the exper- 
iments, Ti is essentially proportional to the DBR's resid- 
ual transmission, and can be taken as constant. Wi^ s 
is the scattering probability from the initial \i) to the 
final |s) polariton state per unit of time. T s describes 
the transmission of the final polariton states to the fi- 
nal photon state outside the sample, which is detected. 
Since the final polariton state leaves the sample due to 
its photonic part, T s is essentially the projection of the 
polariton state on the final photon continuum state, and 
is thus proportional to the mediating polaritons photon 
strength, given by |^4iy/| 2 - 

We will center our attention on Wi-> s , which actually 
describes the inelastic scattering process itself. The tran- 
sition probability is given by Fermi's golden rule as 



Wi- 



2tt 

Y 



\(i\H'\s)\ 2 p(co), 



(9) 



where the initial polariton state is \i) = p\, |0), the fi- 
nal state is given by \s) = p\b^ 0), where W creates a 
LO phonon, and H' is the interaction between polaritons 
and LO phonons (e.g. Frohlich-interaction) through the 
excitonic part of the former. And p{u>) stands for the 
density of final polariton states and is given by 



(10) 



Since the Hilbert space defined by ^ and (J7J) is infi- 
nite, it is cumbersome to work with the eigenstates and 
eigenenergies. However, we do not need them to ob- 
tain the intensity of the Raman scattering process. It 
is enough to use retarded Green's functions, that involve 
the photon and exciton operators (/?/)• These functions 
Gji(ui) — {{Pj\p[))u can be obtained solving a system of 
equations derived from the equations of motion for the 
Green's function: 



(11) 



The polariton-phonon interaction H' is proportional 
to the exciton part of the scattered polariton. Thus, for 
the case of a system where only two exciton states inter- 
act strongly with the MC mode, the Raman intensity of 
Eq. © results proportional to 



I(lo) cx \A y f\ 2 \ A ve i + aA ue2 \ 2 p{uo) , 



(12) 



where v is the label of the polariton eigenstate such as 
E v = lo. 

Here we are neglecting the contribution of the electron- 
hole continuum to H' , and a represents the ratio of ma- 
trix elements of the exciton-LO phonon interaction be- 
tween 2s and Is excitons. Note that if lifetime effects 
are neglected, i.e. if all states but the three polariton 
branches are left aside, p„(lo) is reduced to 5 (to — E v ), 
for v = LP, MP, and UP, and the simple three-mode 
model for the RRS intensity is recovered^ 

Using the Lehman representation of the Green's 
function 33 it can be shown that 



1 



[G 3 i(u 



2tt_' 
A vj A v ip{u) , 



Gjiiu-iQ-)] 



(13) 



and taking into account that Ylj \^j\ 2 — 1) an( l f na f 
p(u>) — J2j Pjji^i it can be followed by replacing these 
expressions in Eq. (TT2"|) , that the desired equation for the 
Raman intensity is given by 



I(u>) oc 



Pff[Pel,el + H 2 |Oe2,e2 + 2Re(a; p e 2,el )] 

E,- Pjj 



(14) 



In the denominator we neglect the contribution of the 
continuum states, because it is negligible near the res- 
onance condition for the outgoing polariton. Therefore 
the sum over j includes only the densities of excitons and 
photons. In particular if only two excitons are involved in 
the polariton system, we need to calculate only the diag- 
onal Green's functions G//(w), G e i je i(o;), and G e 2, e 2(w). 
And also the crossed Green's function G e 2 je i(w). Using 
the equations of motion we obtain: 



G e2 ,e2(w) 



1 



1 



S' 



L0 2 



1 

V? 



G e 2,el(w) = G e l,e2(w) 

V!V 2 



LO\LOiLOj 



V 2 Cj2 — Vz&i 



(15) 



where 



ujf — lo - Ef - Sf , Qi=LO - E et - S ei , (16) 



and 



lo + i0+ - e 



S ei {uj) = ^ 

1 



l^l 2 

LO + i0+ - € 



III 



lo + i0+ - e k 



(17) 
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For the first two of these sums, we assume that the results 
are imaginary constants that we take as parameters: 

S/H = -% , Sj(u) = -ijej ■ (18) 

This is the result expected for constant density of states 
and matrix elements. These two approximations are 
clearly not valid for the last sum. The last sum in 
Eq (fT7|) , corresponds to the excitons that do not interact 
strongly with the MC photons, and to the electron-hole 
continuum. This continuum begins at the energy of the 
gap and corresponds to vertical transitions in which the 
light promotes a valence electron with 2D wave vector k\\ 
to the conduction band with the same wave vector. In the 
effective-mass approximation, the energy is quadratic 
with fen and this leads to a constant density of states 
beginning at the gap. As mentioned in the previous sec- 
tion, Vk ~ {k c \Po\k v ). Approximating the wave functions 
as plane waves and averaging over all directions of k\\ one 
has (|Vfe| 2 ) ~ fcjj/2, proportional to the energy Us- 
ing these assumptions and separating S'Aw) in real [r(u)] 
and imaginary parts we can write: 

S' f (u) = r(w) -<»(«), r(w) = - [ def^L- ,(19) 
1 7T J (u - e) 

y(w) = A(u-Exc)&(u-Exc), (20) 

where Ex c is the bottom of the electron-hole continuum, 
A is a dimensionless parameter that accounts for the cou- 
pling magnitude of the interaction between the MC pho- 
ton and the continuum of excitations, and Q(u>) stands 
for the Heaviside function (step function). The real part 
r(uj) can be absorbed in a renormalization of the photon 
energy and is unimportant in what follows. The imag- 
inary part y(to) is a correction to the photon width for 
energies above the bottom of the continuum. From the 
formalism outlined above, this is expected to lead to a 
decrease in the RRS intensity when the energy is above 
the bottom of the continuum. The expression for y(u) 
accounts so far only for effects of the exciton-continuum, 
and leads therefore to a discontinuous first derivative of 
the calculated RRS profile. This discontinuity is over- 
passed, if the remaining bound exciton states are con- 
sidered. To account for this infinite number of bound 
states, which become very dense below the bottom of 
the continuum, the lower part of y(u>) is replaced by a 
parabolic function that matches (uj — Exc) with a contin- 
uous derivative, starting from zero at the energy (E' x ). 
At the energy E' Xa , which corresponds to the bottom of 
the electronic excitations in H c , as described in the pre- 
vious subsections: 

{0 if uo < E' Xc 

A ^^£l ME' Xc <u<2E Xc -E' Xc . 
A{lo- Exc) if v > 2E Xc - E' x 

(21) 

As it will be show in the next section, this refine- 
ment in the description of the continuum electron-hole 



electronic excitations and weakly bound excitons does 
not affect the Raman profile in resonance with the LP 
and MP branches, and only affects the resonant profile 
on the UP branch for energies above E' x , leading to 
smoother curves. 

The Raman intensity I{u>) for a fixed detuning is given 
by Eq. (H]), together with Eqs. CCH, ([TSJ)— dHJ), (J3TJ) and 
its dependence with detuning (or equivalently with the 
spot position) enters through the corresponding varia- 
tion of the bare MC photon {EA and exciton (E ei ) ener- 
gies. In the next section we will discuss this model, and 
compare the results to the above presented experimental 
data. 



V. ANALYSIS AND DISCUSSION 

The role of the cavity polaritons as intermediate states 
in the scattering process is fully displayed in the reso- 
nant enhancement of the Raman intensity as function 
of exciton-photon detuning, as described in Sec. IIIII and 
shown in Fig. [2] Here the intensity of the scattered reso- 
nant Raman signal tuned to a polariton mode (the MP in 
this case) is followed continuously as this mode changes 
its character from very excitonic (ls-type) to photonic, 
and again very excitonic (2s-type). 

In a simple picture, since the interaction of the po- 
lariton with the LO-phonons is produced through the 
excitonic part of the polariton, the inelastic scattering 
process is favored when the polariton excitonic part is 
large. On the other hand the Stokes shifted polariton 
is only detected when the polariton is transmitted out- 
side of the sample. Since this coupling to the photon- 
continuum is provided through the photonic part of the 
polariton, the Raman process is favored when this weight 
is also large. Thus in view of these two considerations, 
the process is optimized when a compromise between 
both situations is found, and this is reached at zero de- 
tuning, i.e. at the crossing of the bare exciton and MC 
photon, or what is equivalent at the energies where the 
differences Emp — Elp or Ejjp — Emp have a minimum. 
Thus, within this simplified picture, at these energies the 
RRS efficiency is maximized, and the scattered inten- 
sity decreases towards the pure excitonic and photonic 
limits. As we will see, the inclusion of photon and ex- 
citon lifetime broadening, i.e. the coupling of the MC 
photon to the exciton continuum described in Sec. HV| 
introduces important modifications to this correct but 
simplified view. 

For comparison, the results of calculating the RRS in- 
tensity using the simple model given by Eq. (|12p . con- 
sidering only three branches (with no lifetime broaden- 
ing) is shown in Fig. (dashed lines). Two maxima 
are displayed corresponding to the energy of the two 
anticrossings (compare to Fig [2^) . This can be under- 
stood when the corresponding strength^ components of 
the MP are analyzed (plotted in Fig []J>) . The first an- 
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ticrossing (~ 2.2 mm) happens at Ef — E% s = where 
\Af\ 2 ~ | ^4 ls | 2 since the 2s strength at this position is 
very small. The second anticrossing (~ 3.4 mm) hap- 
pens shortly before Ef — E 2s = 0, where the total exci- 
ton part equals the photonic part. As explained above 
at these points the compromise is found, and a maxi- 
mal scattering intensity is expected within this model. 
As one can see in Fig this model qualitatively re- 
produces the general behavior well. However as pointed 
out befor o 18 ' 19 many features are clearly not reproduced 
quantitatively. The first experimental maximum appears 
shifted (~ 2 meV) to higher energies, and something sim- 
ilar happens to the second peak, which appears shifted 
to lower energies (~ 1 meV) and less intense (shoulder 
like), as indicated with arrows in Fig. [2fc. These obser- 
vations can be assigned to lifetime or damping effects^ 
and point to the need of the polariton mediated Raman 
theory described in the previous section (Sec. HVBp that 
includes rigorously the lifetime broadening into the po- 
lariton states. In what follows, our experimental results 
will be compared to this theoretical model based on a 
Green's function approach. 

The Raman intensity, is given by Eq. (fT4")> . together 
with Eqs. flgp. (fT5 ) -([T8" ]) . Parameters Ef, E ls , and E 2s 
as function of the spot position x, as well as 2 V\ s — 
Q,f-i a and 2 Vis = Qf-2s, were derived from the fit of 
the experimental dispersion shown in Fig [T£i, using the 
result of the three-coupled-mode model as a seed. We 
estimated the bottom of the exciton continuum energy 
at Exc — 1.670 meV. Parameter A is taken as 0.031, de- 
rived from the fit to Raman experiments in resonance 
with the UP on a similar II- VI cavity in the very strong 
coupling regime^ The line broadening of the non in- 
teracting cavity-photon is estimated by standard realis- 
tic reflectivity calculation o 36 ' 37 to be 7/ ~ lmeV. The 
exciton line broadening 7 el , 7 e 2, and the 2s-ls ratio 
of the exciton-LO phonon interaction matrix elements 
given by a (assumed as a real number), together with 
the proportionality factor, are left as fitting parameters. 
Both types of experiments described in Sec. IIIII (Fig 
and Fig [4]) are fitted simultaneously. Summarizing in 
a simple interpretation, the first experiment for vary- 
ing spot position changes the photon-exciton character 
of the mediating polariton. While for the second exper- 
iment, the resonance scan at a fixed spot position (fixed 
detuning), accounts for the total polariton line broaden- 
ing for that particular exciton-photon mixture. Best fits 
were obtained for 7 el = 1.2 me V and 7 e 2 = 0.9 meV, and 
a = —0.45. The resulting curves are shown for the RRS 
in Fig. [2J; (full line) and for the resonant scans respec- 
tively in Fig. 2] (full lines). To fit the latter, only the 
upper half part of the scan has been used for the (b) and 
(c) scans, while for (a) the whole scan has been used. As 
can be observed the agreement with the measured data is 
very good, and provides a quantitative description. The 
obtained values for the parameters are relatively robust 
whether a, j e i or j e 2 are changed. 

Since the Is- and 2s-exciton weights for the three po- 




x - spot position [mm] 

FIG. 5: (Color online) Contour plot of the total calculated 
outgoing Raman intensity. The LP and UP mode positions 
are indicated with dashed lines, and the MP outgoing reso- 
nant condition is plotted with a full curve. The Raman scans 
of Fig [4] are indicated with vertical segments and labeled (a), 
(b), and (c) respectively. The intensity is given in logarithmic 
scale. 

sitions [scans (a),(c) and (b)] are different (see FigJTJj), 
small changes in these parameters modify the width of 
the calculated lorentzian-type curves obtained for a fixed 
detuning, changing consequently also the shape of the 
RRS scan. The simultaneous fit of both independent ex- 
periments determines unambiguously the way the values 
that these parameters can adopt. 

A point that might seem odd at first sight, is the fact 
that for the calculations 7 e i results slightly larger than 
7e2- Nevertheless, this result is to be expected if the 
origin of the broadening is disorder with a correlation 
length smaller or of the order of the Is exciton's sizei 29 i 30 
the inhomogeneous broadening is substantially reduced 
by averaging over a volume larger than the length scale of 
the disorder potential. In fact, this is the reason why the 
disorder does not affect the broadening of the photons, 
the wave length of which is much larger than the disorder 
correlation length^ 

The total calculated Raman outgoing intensity is pre- 
sented in Fig. [5] as a contour plot as function of the spot 
position (x) and the corresponding Stokes-shifted energy. 
The mode positions are indicated with dashed lines (LP 
and UP), and the MP outgoing resonant condition is 
plotted with a thick curve. In this figure the Raman 
intensity is shown in a logarithmic scale. The energies 
and x-positions of the Raman scans (shown in Fig. 3|) 
are also indicated with vertical segments [marked as (a), 
(b), and (c) respectively], and the corresponding exact 
resonant position with a full circle. As can be seen in 
this figure the Raman intensity decreases slightly asym- 
metrically when moving away from the exact MP Raman 
outgoing resonant condition when the Stokes-shifted en- 
ergy is varied near the position of the first anticrossing, 
where the LP and MP modes are closer. It has to be 
pointed out that the calculated asymmetry points in the 
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4 1.62 1.63 1.64 1.65 1.66 1.67 

Stokes-shifted energy [eV] 

(b) (a) (c) 




FIG. 6: (Color online) Three-dimensional plot of the calcu- 
lated total outgoing Raman intensity profile as function of 
spot positions (x) and the corresponding Stokes-shifted en- 
ergy. The full curves (LP, MP, and UP) projected on the 
xy-plane of the plot indicate the polariton dispersion of the 
polariton branches. The highlighted resonance profiles on the 
MP, marked as (a), (b), and (c) correspond respectively to 
the fitted curves in Fig. [3] The bottom panel shows a rotated 
detail of the UP region of the upper panel. 



right direction, but is unfortunately not enough to ac- 
count for the stronger one observed experimentally for 
Raman profiles (a) and (c). Something similar happens 
at the second anticrossing between the MP and UP, but 
in this case, due to the extremely low resonant outgo- 
ing Raman intensity on the UP, the anticrossing is not 
observable and is "washed out" . 

The same calculated total outgoing Raman intensity 
profile is presented in Fig. [6] as a three-dimensional plot. 
The mode positions (LP, MP) derived by projecting the 
corresponding resonance maxima on the xy-plane, are 
shown with dotted curves. For the case of the UP mode, 
as previously mentioned, it was not possible to determine 
the precise position of the RRS maximum of the UP at 
the second anticrossing. Thus, for this case, the curve 
shown corresponds to the one derived from the three- 
coupled- mode model (see Fig.QJi), which differs very lit- 
tle from the one derived by projecting the corresponding 
UP resonance maxima (were the UP resonance maxima 
are clearly distinguishable) . The bottom panel of Fig. [5] 



shows a rotated detail of the UP region. The UP reso- 
nant profiles are presented highlighted. Here the upper 
Raman resonant peak is barely seen. For small x, the UP 
appears as a tiny peak that increases slightly its intensity 
when approximating the second anticrossing, and merg- 
ing shoulder like to the resonant profile coming from the 
MP, and increasing its energy after the anticrossing, de- 
creasing its intensity. As can be derived from Fig. [T^, for 
small x the UP has a strong excitonic character, which 
is conserved shortly before the second anticrossing. The 
small photonic component of the UP at these energies 
is very small (see Fig. [5]), and thus the coupling of the 
Stokes-shifted polaritons to the continuum of photonic 
states at the exterior of the cavity is very weak. When the 
position approaches the second anticrossing the photon- 
exciton mixture of the UP becomes larger, and therefore 
the Raman intensity increases, reaching a broad maxi- 
mum slightly after this anticrossing. For increasing ener- 
gies, the coupling of the electron-hole continuum to this 
mode becomes larger and consequently its life time is 
reduced, the RRS intensity decreases, and the Raman 
profile becomes significantly broader. 

As can be seen from Fig. [51 the scattered intensity for 
the outgoing resonant Raman condition with the UP is 
more than one order of magnitude smaller as compared to 
the situation in resonance with the MP. This extremely 
small light intensity scattered by the upper branch near 
the second anticrossing, together with the important in- 
tensity loss due to the interaction of the polariton with 
the exciton continuum (as explained above) are the most 
probable reasons that hindered the experimental obser- 
vation of this mode. For completeness, the fitted Raman 
intensity profiles of Fig. 3] are shown on the same plot, 
also highlighted by thicker lines and respectively indi- 
cated as (a), (b), and (c). 

In our theory, we have not considered the possibility 
of polaritons of a higher branch to be scattered to lower 
polariton states. Since the LP is the lowest excited state 
of the system, the inclusion of this effect would cause an 
important increase of the LP resonant outgoing Raman 
intensity, with respect to the resonant scattered intensity 
by the higher order polariton branches (MP, UP). This 
kind of effect has been reported by Tribe et al^ 1 for a 
III-V type cavity with build-in electric field. The signifi- 
cantly larger contribution to the total scattered intensity 
coming from the LP, may explain the asymmetry of the 
observed Raman scans (see Fig. HJ, and would be partic- 
ularly important at those situations were the LP and the 
MP are nearest, as is the case in for the Raman scans 
labeled (a) and (c) in the Figs. 0] to [H i.e. near the first 
anticrossing. 

A study of this processes, as well as detailed calcu- 
lations of the different contributions to the exciton line 
width (supposed here as simple constants), and the care- 
ful determination of the exciton-LO phonon matrix ele- 
ments (here accounted by the constant factor a), might 
lead to refinements and improvements of the theory. 
However, we believe that the formal framework of how 
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the main physical ingredients enter the polariton medi- 
ated Raman process, is well described by the above pre- 
sented theory. 

VI. CONCLUSIONS 

Summarizing, we have reported two types of indepen- 
dent Raman experiments in near resonance with the mid- 
dle cavity-polariton branch of a II- VI based semiconduc- 
tor optical cavity. The first experiment shows the behav- 
ior of the Raman intensity at exact outgoing resonance 
with the cavity-polariton, for varying photon-exciton de- 
tuning. And the second experiment shows the evolution 
for a fixed photon-exciton detuning, when the Stokes- 
shifted photons are tuned on and off resonance by varying 
the excitation energy. Both independent experiments are 
accounted simultaneously by the here developed cavity- 
polariton mediated Raman theory, that rigorously in- 
cludes photon and exciton damping effects. The pre- 
sented theory describes quantitatively well the experi- 
mental results, a fact that remarks the importance of the 
exciton and cavity-photon lifetime effects in the process 
of inelastic scattering of light due to QW optical phonons. 



The shift of the Raman intensity maxima towards the sit- 
uation where the polariton photonic strength is larger is 
one of the important features that are predicted by this 
model, and reasonable values for the damping constants 
are derived. 

The theory in its present form is not able to explain the 
asymmetry of the observed Raman scans (see Fig. |4]). It 
might be possible that inelastic mechanisms of decay of 
the middle polariton in the lower one, not included here, 
might lead to an improvement of the comparison with 
experiment. Nevertheless, we believe that the presented 
Raman scattering theory based on a Green's function ap- 
proach, allows a consistent introduction of the damping 
mechanisms, and provides an important step forward in 
understanding this fundamental and interesting aspect of 
the inelastic light-matter interaction processes in optical 
microcavities. 
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